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Abstract 

We discuss whether quark, charged lepton and neutrino masses and mixing angles may 
be related by an extended flavour and family symmetry group. We show that current 
measurements of all fermion masses and mixing angles are consistent with a combination 
(2/)! of an underlying SU(3) family symmetry together with a GUT symmetry such as 50(10). 

In this the near bi-maximal mixing observed in the neutrino sector is directly related to the 
small mixing observed in the quark sector, the difference between quark and lepton mixing 
angles being due to the see-saw mechanism. Using this connection we make a detailed 
prediction for the lepton mixing angles determining neutrino oscillation phenomena. 

> : 

1 Introduction 

<N ; 

While the measurement of quark and lepton masses and mixing angles continues to improve, 
their theoretical understanding remains elusive. In the Standard Model the quark masses and 
mixing angles and charged lepton masses are determined by the Yukawa couplings, parameters 
which must be specified when defining the Standard Model. The renormalisable couplings of 
the Standard Model do not allow neutrino masses although they can be introduced through 
Ph! the addition of higher dimensional operators, presumably originating in physics Beyond the 
Standard Model. 
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There are a few promising suggestions for structure Beyond the Standard Model which ad- 
dress the fermion mass problem. For simple symmetry breaking schemes Grand Unification can 
relate quark and lepton masses. The most promising of such relations is the equality m& = m T , a 
result which applies at the GUT scale. Radiative corrections are dominated by QCD interactions 
which increase the bottom quark mass at low energy scales, giving a prediction for the bottom 
quark mass in good agreement with experiment. A fairly straightforward generalisation of such 
GUT relations, which invokes a new family symmetry, also provides good relations between the 
down quarks and charged leptons of the first two generations [0. 

However the recent measurement of neutrino masses and mixing angles has shed doubt on 
the validity of simple relations between quark and lepton masses and mixing angles. The most 
obvious difference is the fact that neutrinos have masses much smaller than those of the quarks or 
charged leptons. In the context of Grand Unification the "see-saw" mechanism |3|, f|] provides 
an elegant mechanism to explain this difference. In this scheme the Standard Model spectrum 
is extended to include the right-handed (SU (2) singlet) partners of the SM neutrinos. As these 
states can acquire SU(3) x SU(2) x £7(1) invariant Majorana masses, naturalness implies they 
should acquire mass at the scale (the GUT scale?) at which the theory Beyond the Standard 
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Model breaks to the Standard Model. In addition one expects Dirac masses between the left- 
and right- handed neutrino states. These are generated by electroweak symmetry breaking in 
a similar way to that of the charged leptons and quarks and so it is reasonable to assume they 
will be of comparable magnitude. Diagnosing the matrix of Dirac and Majorana masses gives 
the mass matrix for the three light states, dominantly composed of the doublet neutrinos, of the 
form 



Here M v is the 3x3 matrix of Majorana masses for the three generations of right-handed 
neutrinos and M£ is the 3x3 matrix of Dirac neutrino masses. Given that the natural size of 
M£ is much less than that of M u , it is clear this offers an elegant explanation for the smallness 
of neutrino masses. 

There remains the question why lepton mixing angles should be so different from quark 
mixing angles. This is the question we address here in the light of new information about the 
quark mass matrix that has recently been obtained from an analysis of the wealth of b-quark 
data and the information from neutrino oscillations. Our approach is to look for the largest 
symmetry consistent with the fermion masses concentrating on the possibility that quark and 
lepton masses are related. We show that the data in the quark sector data is consistent with 
an underlying (spontaneously broken) non-Abelian family symmetry together with an extension 
of the Standard Model flavour symmetry to include an SU(2)r flavour symmetry acting on 
the right-handed components. The charged lepton masses and neutrino masses and mixings 
are consistent with this symmetry provided the flavour symmetry is extended, for example to 
SO (10), to generate GUT relations between lepton and quarks. This structure can lead quite 
naturally to large lepton mixing angles. Indeed, due to the see-saw mechanism, lepton mixing 
angles are large because quark mixing angles are small || Q ! 

2 The quark mass matrix 

Our starting point is the form of the quark mass matrix. Since experiment is only able to de- 
termine the mass eigenvalues and the CKM mixing matrix, it is not possible unambiguously to 
determine the quark mass matrix because the full form of the left- handed and right-handed 
rotation matrices needed to diagonalise the quark masses is not known. However the very rea- 
sonable assumption that the smallness of the mixing angles is due to the smallness of the mixing 
angles in the up- and down- sectors separately, allows one to determine the mass matrix elements 
on and above the diagonal to good precision for the down quarks and to lesser precision for the 
up quarks. In reference 0, under this assumption, the quark mass matrices were determined 
using the most recent experimental data. This gave the form 
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In this the parameters of the up quark mass matrix are given by e = 0.05, V ~ e 1 ^, where in the 
phase convention used here « 90° is the 'standard' (i.e.PDG convention) CP violating phase S 
|7j while a' and d are very weakly constrained. The parameters of the down quark mass matrix 
are much better determined with 

e = 0.15 ±0.01 b = 1.5 ±0.1 a =1.31 ±0.14 
|c| = 0.4 ±0.05 ip = Arg(c) = -24° ± 3° or 

|c| = 1.28±0.05 ^ = 60° ± 5° (4) 

In both M u and M d the entries marked with a question mark are only weakly constrained. In 
what follows we shall consider the case that the mass matrices are symmetric or antisymmetric 
or some combination of the two so that the elements below the diagonal are also determined. 
This is in the spirit of looking for the maximal symmetry consistent with the masses and is 
motivated by certain GUTs, notably 5O(10). There is one strong piece of experimental evidence 
for such a symmetric structure, namely the success of the Gatto, Sartori, Tonin relation || 
between the Cabibbo angle and the quark masses of the first two generations 




where a is the CP violating phase entering the Jarlskog invariant |7|, |9| . This relation only applies 
if the (1,2) and (2, 1) matrix elements are equal in magnitude. The solution of eqs@ and (||]) 
has this structure. Of course there is no direct indication that the same reflection symmetry 
applies to the remaining matrix elements but it is the most natural generalisation of this result. 
In any case, as stressed above, we think it of interest to ask whether the most symmetric form 
for the quark mass matrices can be simply related to viable forms for the lepton mass matrices. 



3 Quark family and flavour symmetries 

In this Section we will discuss whether the quark mass matrices are consistent with a larger 
set of symmetries than those of the Standard Model^J These may be flavour symmetries acting 
in the same way on each of the three families; for example simple Grand Unified theories are 
flavour symmetries. Alternatively they may also be family symmetries distinguishing between 
families. In our opinion the structure of the quark mass matrices strongly suggests an underlying 
(spontaneously) broken family symmetry. In such a scheme, in leading order, only the (3, 3) 
elements are allowed. In terms of the parameter, e, characterising the symmetry breaking 
the remaining matrix elements are filled in at some order e n which is also determined by the 
symmetry, thus generating the hierarchical structure of fermion masses. 

The candidate family symmetry group, namely those symmetries which commutes with the 
gauge interactions of the Standard Model, is quite large. For the quarks it is Z7"(3) 3 , corresponding 
to a separate Z7 (3) factor for the left-handed doublets and each of the right-handed SU{2) singlet 
fields respectively. To generate quark mass matrices with a left- right- symmetry we restrict the 

3 Here we shall only consider supersymmetric theories. 
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choice to one with the same transformation for the left-handed and the charge conjugate right- 
handed quarks, this symmetry is reduced to a diagonal £7(3) or one of its subgroups^. 

Abelian family symmetries have been widely explored. They are capable of explaining the 
hierarchical structure of the quark masses and generate quite acceptable forms for the quark 
matrices [K], j5|. The main shortcoming of these schemes is that they do not give a reason for 



the smallness of V c b- This follows because if one assigns Abelian charges to generate the mass 
hierarchy between m s and m^, it leads to the prediction for V& — ^/mjrrib which is too large 
HT0| , |Tl|j . Although this relation is quite sensitive to corrections of 0(1) due to unknown Yukawa 
couplings, it is an unsatisfactory feature of an Abelian symmetry. To obtain the correct form for 
V c b it is necessary for the ratio of the (2, 2) and (2, 3) matrix elements to be close to 1; a ratio 
of 1 implies V c b — m s /mb, in good agreement with experiment. In contrast to the Abelian case 
a non- Abelian family symmetry can relate the magnitude of the couplings in different positions 
in the mass matrix and thus offers the possibility of explaining this near equality. 



3.1 Non- Abelian family symmetry 

The largest non-Abelian symmetry which treats the left- and charge conjugate right- handed 
components in the same way is SU(3). In [H] a specific SU(3) theory has been constructed which 
relates the (2, 2) and (2, 3) matrix elements. In it the left-handed quark doublets and the up and 
down charge conjugate quarks are assigned to SU(3) triplets, ipi, u\ and d\ respectively. The 
critical part of the theory lies in the pattern of symmetry breaking which leads to the desired 
mass matrix. This is achieved through the antitriplet scalar fields 023; 03 an d triplet fields 2 3j, 
4> 3i which acquire vacuum expectation values (VEVs) 03 = ( M ) , 0^ = ( a 2 a 2 ) , 
03 = ( M ) and 23 = ( a 2 —a 2 ) . In || it was shown how such vacuum alignment 
can be achieved in a supersymmetric theory, with M > a 2 . 

With this pattern of symmetry breaking it is easy to generate the required form of the mass 
matrices via the superpotential 

PYukawa = (#/4^+#/>230^ (6) 

where H 12 are the Higgs doublets of the MSSM, which are SU(3) singlets. If ip c — u 1i dt 
then a = 2, 1 respectively, giving rise to the up and down quark masses in a supersymmetric 
theory. In eq(|J) we have suppressed the couplings of 0(1) associated with each operator. The 
first two operators shown are the leading (dimension 6) ones consistent with a simple family 
symmetryQ[0. The third and fourth operators are only generated via higher dimension operators 
involving additional powers of fields 023, 03. Replacing these fields by their vevs leads to the 
suppression factors e 2 and e 6 shown, where e = a 2 /M. This gives rise to the mass matrices of 

4 It is possible there is an even larger symmetry with separate gauge group factors operating on left-and right- 
handed components and related by a discrete reflection symmetry. If one requires the left- right- symmetry is 
preserved on spontaneous symmetry breaking, the phenomenological implications of this scheme are the same as 
the case considered here. 

5 In a specific realisation of the scheme it is necessary to introduce further Abelian family symmetries to ensure 
this is the most general form. A specific realisation is given in reference [[| . 
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the form 



M 



M: 



3,3 




(7) 



where the expansion parameter e may differ for the up and down sectors and we have explicitly 
included the 0(1) coefficients A, A', A" that arise due to the 0(1) couplings in eq(|) which are 
not related by the non-Abelian symmetry. The effect of the 577(3) family symmetry may be 
readily seen relating the (1,2) and (1,3) matrix elements as well as the (2,2) and (2,3) matrix 
elements. Of course there are corrections to these equalities coming at higher order in e from 
operators of higher dimension. We shall consider their effect in detail later. 

The form of eq(0) gives good agreement with the measured quark mass matrices and can 
also be extended to give a viable description of charged lepton and neutrino masses and mixing 
angles 0. Here we consider whether the 0(1) coefficients can be related by further symmetries 
and whether it is possible to determine the lepton masses and mixing angles in terms of the 
quark masses and mixing angles. 



3.2 SU(2) R and Grand Unification 

In implementing the SU(3) family symmetry the family structure has been strongly constrained 
to generate left-right symmetry by assigning the left- and the charge conjugate of the right- 
handed components of a given quark to the same representation. Since, in addition, the family 
symmetry must respect the SU(2)l gauge symmetry the up and down left-handed components 
of a given quark family doublet must also carry the same family charge. As a result of these 
two constraints one immediately sees that left- and charge conjugate right- handed components 
of the up and down quark members of a given family have the same family charge. This means 
that the up and down quark mass matrices have the same form as in eqs(|^) and (|3]), although 
the expansion parameters may be different and the operator coefficients of 0(1) may differ. 

If SU(2)r is also an exact symmetry of the theory even the expansion parameters and 
operator coefficients are equal and the mass matrices are identical. This is clearly not acceptable 
and, if there is an underlying SU(2)r symmetry, it must be spontaneously broken so that the 
equality of the mass matrices will be lost through soft symmetry breaking terms. These enter 
through the expansion parameter e which is determined by 9/M where 9 is the field spontaneously 
breaking the symmetry and M is the messenger mass of the state responsible for communicating 
the symmetry breaking and generating the higher dimension operators. Due to SU(2)r breaking 
the messenger mass may be different for the up and down quark sectors and hence the expansion 
parameters may differ. It is also possible that the family symmetry breaking field, 9, is not a 
singlet under SU (2) R and its vev breaks SU(2) R , again leading to a different expansion parameter 
for the up and the down sectors. However there will still be some measurable effects of an 
underlying structure SU (2)r in the quark masses because the Yukawa couplings of the field 9 to 
the quarks and the messenger sector, which are responsible for the 0(1) coefficients, will respect 
the symmetry. How large are the corrections to this result? They occur in radiative order but, 
in a supersymmetric theory, supersymmetry guarantees that superpotential is not renormalised. 
In this case the dominant radiative corrections to the Yukawa couplings come from the D-terms 
controlling wavefunction renormalisation. Thus we may expect that fields falling into SU(2)r 
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doublets will have normalisation which differ only in higher order, 0(\6\ 2 / M 2 ). The resulting 
corrections to the 0(1) coefficients will also be at this order. 

What are the phenomenological implications of an underlying SU(2)r? If only the leading 
operator contributions to eq(0) are included the coefficients a, b and c in eq(§) will equal a', b' 
and d in eq(|2|) respectively^]. We wish to explore whether such a broken symmetry expansion is 
consistent with the quark and lepton masses and mixings. Unfortunately the tests of the SU(2)r 
relations in the quark sector are somewhat imprecise at present because the contribution to the 
mixing angles from the up sector is small and sensitive to small changes in the down sector. As 
a result a 1 and d are only poorly determined. Given that in eqs(^) and ([J) the (2, 2) matrix 
elements have been used to define the expansion parameters, the remaining test is to compare b 
and V . It is easy to eliminate the dependence on the expansion parameters by combining (1, 1), 
(2,2) and (1,3) elements. Using eqs(|2|) and (^) this gives 

b' m s lm u m t , , 

= —4/-^ (8) 
b m c y mdrrib 

where this relation applies at the unification scale where SU (2) R is a good symmetry. This gives 

^- = 0.5 ±0.3 (9) 


Given that the expansion parameter in the down quark sector is quite large, e = 0.15, and 
that higher order corrections may be significant, the result of eq(||) may be consistent with an 
underlying SU(2) R relation prediction between b' and b, with a significant correction at 0(6/ M) 
coming from such higher dimension operators. We shall return to a discussion of this possibility 
in Section |j. 

Of course it is likely that SU(2) R is part of an underlying GUT such as 50(10). In this 
case the question whether SU(2) R is an approximate symmetry of the couplings depends on 
the pattern of symmetry breaking. We shall also discuss a second possibility for reconciling the 
GUT prediction with eq(||) i n which the couplings themselves feel SU(2) R breaking. 

3.3 GUTs and charged lepton masses 

We have seen that the quark mass matrices are consistent with the choice that all the states of the 
(left-handed) components of a given family have the same transformation properties under the 
family symmetry. This is suggestive of a larger underlying GUT symmetry. The GUT 5*0(10) is 
particularly promising as all the (left-handed) states of a family, plus the charge conjugate of the 
right handed neutrino, fit into a single 16 dimensional representation. If 5*0(10) is an underlying 
symmetry of the theory any family symmetry must commute with it implying all the charges of 
a given family must be the same, consistent with the form of the mass matrix discussed above. 
An associated advantage of such a family symmetry is that the mixed anomalies of the Standard 
Model gauge group with the family symmetry will automatically cancel because of the structure 
of the underlying 5O(10). It also contains SU(2) R and can relate the up to the down sectors as 
discussed above. In addition a GUT symmetry can relate quark and lepton mass matrices and 
this is the issue we wish to study here. 

6 In general we expect higher order operators contributing to a given matrix element at 0(9/ M) or above 
relative to the leading term. 
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We first consider the charged leptons. Given the same family properties the form of the 
charged lepton mass matrix will be the same as that of the down quark matrices in eq(|3|) up 
to coefficients of 0(1). If there is an underlying GUT the coefficients too may be related. As 
mentioned above the relation mb = m T at the unification scale is in good agreement with the 
measured masses. Such an equality applies in SU(5) if the Higgs responsible for the third 
generation masses transforms as a 5 of SU(5). In 5*0(10) equality applies if the Higgs belongs 
to a 10 representation, this also gives equality between the top quark and the third generation 
Dirac neutrino mass, something we explore below. 

What about the two lighter generations? In this case we must address the question whether 
the expansion parameters are related. From a phenomenological point of view, note that, af- 
ter taking radiative corrections into account, the relation Det[M d ] = Det[M l ] at the unifica- 
tion scale is also in good agreement with the experimental measurements. From eq(^) we see 
Det[M]/Mz3 = A 2 e 6 and so equality requires that the (1,2) matrix element of magnitude Ae 3 
be the same for the down quarks and the leptons. This is consistent with an underlying broken 
SU(2) R symmetry because the down quarks and leptons are both T R3 = — 1/2 states and both 
can acquire their mass from the same Higgs doublet, H 2 , in a supersymmetric theory. Thus the 
strong SU(2) R breaking, needed to split the T R3 = ±1/2 states, is consistent with this equality. 

3.3.1 Symmetry breaking expansion parameters 

Of course equality of the down quark and charged lepton matrix elements in the (1,2) position 
requires that the expansion parameters be the same in the two sectors. This is consistent with an 
underlying broken SU (2) R symmetry because the down quarks and leptons are both Tr j3 = — 1/2 
states and both can acquire their mass from the same Higgs doublet, if 2, in a supersymmetric 
theory. Thus the strong SU{2) R breaking, needed to split the T R>3 = ±1/2 states, does not 
differentiate between the down quarks and charged leptons. Clearly if the dominant messenger 
sector for family symmetry breaking is in the Higgs sector the expansion parameter in these 
two sectors will be equal. Similarly, in this case, the expansion parameter in the up quark and 
neutrino sectors will be the same, since both acquire Dirac masses from the same Higgs doublet, 
H\. A similar conclusion applies if the family symmetry breaking field, 8, is not a singlet under 
SU{2) R and its vev breaks SU(2) R , again leading to a different expansion parameter for the up 
and the down sectors. 

The other possibility is that there are significant contributions from messengers carrying 
quark and lepton quantum numbers. For the messengers carrying left-handed generation quan- 
tum numbers, SU (2) L requires the messenger masses should be equal so the only way that these 
terms could be consistent with the up and down masses is if the family symmetry breaking 
field, 9, is not a singlet under SU (2) R and its vev breaks SU {2) R . Finally it is possible that the 
messengers carry the quantum numbers of the left-handed charge congugate generations. In this 
case the up and down sectors could, through SU{2) R breaking, have a different masses leading 
to different expansion parameters. If the underlying symmetry breaking pattern is 

50(10) -> SU{A) x SU{2) L x SU{2) R -> SU(4) x SU(2) L x 17(1) (10) 

the down quarks and charged leptons will have the same expansion parameter and so will the 
up quarks and neutrinos. 
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In what follows we will explore the most predictive possibility, namely that there is an 
expansion paramete e which applies to the down quarks and charged leptons, and an expansion 
parameter e which applies to the up quarks and neutrinos. 

3.3.2 Quark/lepton couplings 

If the expansion parameter is indeed the same for the down quarks and leptons the quantities 
Det[M d ] and Det[M l ] occur at the same order as is requires phenomenologically. Exact equality 
requires also the equality of the coefficients determining the (1,2) matrix elements. Just as for 
the (3, 3) elements such equality will follow if the Higgs responsible for this element transforms 
as a 5 of SU(5) or to a 10 dimensional representation of SO (10). However it is not possible 
to have identical charged lepton and down quark mass matrices because, after taking account 
of the radiative correction on going from high to low scale, giving approximately a factor of 3 
increase in the quark masses, the relations m s = m M and rrid = m e are in gross disagreement 
with experiment. As pointed out by Georgi and Jarlskog |] this discrepancy is readily explained 
if there is an underlying GUT through the appearance of Clebsch Gordon factors in the matrix 
element coefficients. In particular if the Higgs responsible for the (2, 2) matrix element should 
belong to a 45 of SU(5) (or 126 of SO (10)) the lepton coupling is a factor —3 times the down 
quark coupling. In this case, taking account of the equality of the determinants, the relations 
for the light generations are modified to give m s = m M /3 and = 3m e . Including the radiative 
corrections needed to determine the masses at laboratory scales, these relations are in excellent 
agreement with the measured masses. If there is also an underlying non-Abelian family symmetry 
relating the matrix elements as in eq (|3.1|) , the (2,2) and (2,3) matrix elements are related and 
so it is necessary that the Higgs field responsible for the (2, 3) matrix element be also due to 
a Higgs the 45 of 577(5). With this the resulting form of the lepton mass matrix coming from 
eq([|) is given by 



Note that in order to realise this scheme there must be a very particular origin to the effective 
Lagrangian given in eq(|6|). The effective operators may be generated through quark or Higgs 
mixing with states carrying different family symmetry quantum numbers. In order to generate 
the factor 3 it is necessary that the operators responsible for these terms be dominantly given by 
Higgs mixing of the 5 Higgs responsible for the (3, 3) element with the 45 generating the factor 
—3 enhanced terms. Note that it is not necessary for the 45 to be an elementary Higgs field, it 
can arise as an effective Higgs, for example from the coupling of a 5 and 24 (in SO (10) the 126 
can be generated by a 10 and 45) as in Figure 1 [1T2| . In these graphs x an d X (x c an d X c ) are 
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1 



a vectorlike pair of chiral supermultiplets with mass M (M'), where x (x c ) has the same SM 
quantum numbers as the left-handed generations if) (left-handed antigenerations if>°). 

Since the adjoint fields are already needed to break SU(5) (or 5*0(10)) this represents a 
simplification of the Higgs sector. Moreover we see that such graphs can also generate a coupling 
to an effective 120 dimensional representation of 5O(10). This coupling would be antisymmetric 
in family space if the 120 were fundamental but need not be so in the case they are given by 
Figure 1 because the intermediate messenger masses associated with the graphs coupling to H w 
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and E 45 are not necessarily the same as those coupling to S 45 and Hi . Writing the vacuum 
expectation value 

< E >= B - L + K T R>3 , (12) 
the relative contribution of these graphs to the down quarks and leptons is respectively 1/3 and 
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(b) 



Figure 1: Froggatt Nielsen supergraphs generating fermion masses. 

-1 for Figure la and (1/3 — k/2) and (—1 — k/2) for Figure lb. Consider the case the second 
graph has lighter messengers and dominates. For k = we obtain the form of eq([TT|). However 
for case k = 2 the contribution to the leptons is +3 times that of the quarks so we obtain 
eqflll]) with +3 rather than —3. Since the lepton mass eigenvalues are insensitive to this sign 
we obtain identical masses for either form. This case corresponds to the coupling of an effective 
120 dimensional representation and has different implications for the structure of the up quarks 
and neutrinos that we explore below. 

The mass matrix of eq([ll|) gives excellent relation between the down quarks and charged 
lepton masses of all three generations. Due to the approximate texture zero in the (1,1) position 
it also implies a contribution to the (1,2) element of the matrix V 1 diagonalising the charged 
lepton mass matrix, given approximately by 



(13) 



In turn this gives a contribution to the (1,2) element of the MNS mixing matrix, V MNS = V l ^V u 
T3|. As we shall discuss this is significant in solar neutrino oscillations if the (now disfavoured) 



SMA solution applies; it is also important in determining . On the other hand, due to the 

smallness of the (1, 3) matrix element, the contribution of the (2, 3) element of V to the MNS 
mixing matrix is small 



V 1 



(14) 



This is the analogue to the relation in the down quark sector V c b ~ m s /mb, which is in good 
agreement with the small value observed. However it is too small to explain the large mixing 
angle observed in atmospheric neutrino oscillation. This discrepancy is at the heart of the 
apparent conflict between quark and lepton masses and we turn now to the discussion how this 
conflict may be resolved. 
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4 Neutrino masses. 



4.1 Dirac mass 

Following the discussion in Section |3| we immediately see that the SU(3) family symmetry 
properties of the left- and right- handed neutrinos must be the same as those of the charged 
leptons. As a result the form of the neutrino Dirac mass matrix between the doublet neutrinos 
and the singlet (right-handed) neutrinos should be the same as eq ([Tl|) , although the expansion 
parameter may be different. Given the success of GUT relations in the charged mass matrix it 
is obviously reasonable to explore the possibility that neutrino masses are similarly related to 
quark masses. As discussed above the expansion parameter for neutrino masses may be equal to 
that for up quark masses. This is consistent with an underlying SU(2) R (contained in SO (10)) 
because the up quarks and right-handed neutrinos are both T 3 R = 1/2 states. In this case the 
resulting Dirac neutrino mass matrix is of the form 

b' v e* d v e> 

M V D /M» = I b' v e> d> v ? <e 2 I (15) 




The underlying GUT symmetry relates the renormalisable couplings of the messenger states. 
For the case the Higgs responsible for the (3,3), (1,2) and (2,1) matrix elements transforms 
as the 10 of 5*0(10) and the Higgs responsible for the (2,2), (2,3) and (3,2) matrix elements 
transforms as the 126 of 50(10) we have a! v = —3a, d' v = —3, b' v = b, and d v = c. In the case just 
discussed where the Higgs responsible for the (2,2), (2,3) and (3,2) matrix elements transform 
as the 120 of 5O(10) (with the second graph of Figure 1 dominating) we have a' v — — : u ~" ! ' ~" 



H _ (6-3«) 7 / 7 1 / 

v (2+3k) ' v u i aliu 



(2+3k)' 



4.2 Majorana mass 

Of course it is necessary to determine the Majorana mass matrix before one can determine the 
effective neutrino mass matrix via the see-saw formula of eq(|I]). Although the family symmetry 
properties of the right-handed neutrinos are related to those of the charged leptons it is not 
possible to use this information unambiguously to determine the structure of the Majorana 
mass matrix. In particular it may not have the same form as is found for the Dirac matrix. The 
reason is twofold. Firstly the Majorana masses are generated via a new AL = 2 lepton number 
violating Higgs sector and it is necessary to specify the family symmetry representation content 
of this sector before the Majorana mass structure is fixed. Secondly the Majorana mass matrix 
involves the coupling of identical fermions and so antisymmetric terms allowed in the Dirac mass 
matrix will not arise in the Majorana matrix. Despite this we can make some general statements 
about the structure. In particular we expect an hierarchical structure for the Majorana mass 
matrix because the underlying family symmetry (SU(3)), is the same as applies to the Dirac 
matrix which leads to a structure ordered by a new expansion parameter Em- If a single AL = 2 
(effective) symmetry breaking field, $ with definite £77(3) family transformation properties, 
dominates there will no possibility of degeneracy between matrix elements. For a large part of 
the parameter space, the form of the Majorana mass matrix is in practice determined. In the 
limit that em << e (which is the case if the messenger sector in the $ sector is heavier than in 



9 



the electroweak breaking sector) the mixing in the neutrino sector is dominated by the mixing 
coming from the Dirac mass. If $ couples dominantly to the (3, 3) position the resulting mass 
matrix can be diagonalised by small rotations of 0(em) giving 

(mi \ 
m 2 (16) 
m 3 J 

It is important to note that this mass eigenstate basis is likely to be very close to that used for 
the Dirac matrix, eqflTo]), the mass eigenstates will be the family symmetry eigenstates up to 
corrections of 0(em)- 

This is the general form we shall use in our phenomenological analysis. It corresponds to one 
of the two most reasonable possibilities. By "reasonable" , we mean that large mixing should not 
be due to a detailed correlation between the Dirac and Majorana mass matrix elements because 
this would require a correlation between the AL = 2 and AL = symmetry breaking sectors and 
there is no obvious symmetry that could achieve this. The other reasonable possibility is that the 
mixing is dominated by the mixing in the Majorana mass matrix (the case em ^> e).This is less 
predictive (and hence less interesting) because, in this case, we cannot relate neutrino mixing 
angles with quark mixing angles. For this reason we concentrate here on the first possibility. 



5 The light neutrino mass matrix 

We are now able to determine the masses and mixing angles of the light neutrinos using eqs flI5| , |l6D 
in eq(Jl|) . The effective Lagrangian associated with the see-saw mass is of the form 

e 6 e 4 1 

C = — + c> T + 0(e 5 )v e ) 2 H 2 + — + a> T + eb' u u e ) 2 H 2 + — (u T + a' v u v + ec' u u e ) 2 H 2 
mi m 2 m 3 / 17 x 



5.1 Near maximal mixing 

It is now straightforward to determine the viable forms of the mass matrix. Consider first 
the heaviest neutrino which should have near maximal mixing to explain atmospheric neutrino 
oscillation. One sees from eq(]l7|) that, surprisingly, this is easy to achieve. 

We first consider the case — >—,—. The heaviest eigenstate is v\ oc d'u u + a'u T + eb'u e . If 

ni2 mi ' m,3 ° x v H> v < v e 

SU(2)n is an underlying symmetry of the couplings (k = in eqfll^)), a! v /d' v ~ a ~ 1.3. In this 
case v\ is very close to a maximal mixed state of muon and tau neutrinos. Note, as remarked 
in Section |3|, it is the smallness of V c b which requires a (2,3) entry of 0(e 2 ) that leads to near 
maximal mixing in the neutrino sector 0! 

It is also straightforward to compute the next lightest neutrino state. It will have mass of 
0(^-j- or ^) whichever is the larger. Let us consider the first possibility first. The eigenstate 
v 2 has the form 

i/ 2 oc (a' u b' u - c' v d' v )(d' v v T - a' u u^) + 0(e)u e (18) 

Except for the case the factor {a! v b' v — c' u d' u ) is anomalously small the second eigenstate will be 
predominately in the u^, v T direction. The same conclusion applies if the second eigenstate is 
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dominated by the third term of eq(|T7]). Thus in both cases the mixing angle V Vetl relevant to 
solar neutrino oscillation will be dominated by the contribution of 0(e) from the charged lepton 
sector given in eq(pT3[). This is of the correct magnitude to fit the SMA solution. The ratio of the 
two heaviest masses is given by M3/M2 = 0(e 2 m2/mi or e 4 m2/m3) and the Majorana masses 
can be adjusted to get the value needed for the SMA solution. For the case — > — , — a similar 
pattern emerges except that the heaviest state is now v x oc b^v^ + c' v u T . SU(2) R symmetry leads 
to the relation b' v jd v ~ b/c ~ 0.5 and so the mixing is large but not maximal. Provided the 
factor {a' u b' v — d u d' u ) is not anomalously small the lighter states are not strongly mixed and the 
dominant contribution to Vjrf comes from the charged lepton sector. From this we see that a 
large atmospheric neutrino mixing angle and the SMA solar neutrino solution is consistent with 
a very symmetrical structure for the mass matrices in which there is an approximate up-down 
SU(2)l x SU(2)r symmetry together with a GUT symmetry relating quarks and leptons. 

This analysis has demonstrated that, through the see-saw mechanism, it is quite natural 
to get near maximal mixing in the neutrino channel for atmospheric neutrino mixing, the large 
lepton mixing being due to the small mixing in the quark sector! However inclusion of the recent 
SNO data now strongly disfavours the SMA solution for the solar mixing angle and favours large 
mixing in this case too. To address this we now turn to the question whether near bi-maximal 
mixing is consistent with an enlarged family and flavour symmetry. 



5.2 Near bi-maximal mixing 

At first sight is seems that our analysis has already ruled out the possibility of near bi-maximal 
mixing. The possible exception to this is if the factor (a' v b' u — d u d' u ) is of 0(e), in which case 
the second state of eq (|l8D has a large v e component. There are two ways that this may happen 
naturally. The first is through SU(2)r breaking in the coefficients which occurs through spon- 
taneous breaking via the vev of eq(|T2|) for the case k ^ 0. For the special case k — 2, which 
generates good charged lepton masses, we have a' v = d' v = 0. In this case higher order terms in 
the expansion in terms of the symmetry breaking parameter can generate (a' u b' u — d u d' u ) of O(e) 
as required. The case k = 2 can arise quite naturally on spontaneous symmetry breaking as 
it corresponds to an enhanced symmetry point at which the effective potential can easily have 
a minimum. The second way (a' v b' u — d v d' v ) may naturally be of O(e) applies even to the case 
k = and again relies on higher order terms in the expansion in terms of the symmetry breaking 
parameter. 

Let us start with a discussion of the sub-leading terms. Consider the case of the SU (3) family 
symmetry which leads to the leading order mass matrix of the form given in eq(|7]). Including 
higher order terms in the symmetry breaking expansion parameter gives a mass matrix of the 
form 



M/M 3i3 = ( s 3 (z + (x + y)e) e 2 

-e 3 (z + (x — y)e) e 2 {mr 



e 3 (z + (x + y)e) e 3 (z + (x — y)e) 



aw + ue) 
us) 



[aw — ue ) 
1 



(19) 



Here z, w and u are real coefficients and x and y complex coefficients of order 1. For the 
case of the down quarks and leptons the expansion parameter is e = e and for up quarks and 
neutrinos e = e (c/.the discussion of Section |3.3.1|) . The higher order terms are necessary to fit 
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the data for M d of eq(g). Following the discussion of Section |3.3.2| we assume that the coefficient 
a arises from an effective 120 or 126 SO (10) representation in order to generate the correct 
lepton mass spectrum. Thus a is determined by the value of k in eq(|T2"D and there are two 
possibilities 

k = 0, a\ = —3, a u = —3 (20) 
K = 2, ai = +3, a u = (21) 

We assume all other coefficients come from a 10 of 5*0(10) so they are the same as for the 
quarks. 

We may determine the neutrino effective mass matrix using the see-saw formula, eq(|l]). We 
focus on the case that the dominant contribution to the heaviest of the light neutrino states is 
due to the exchange of the right handed neutrino with mass mi which couples to the first row 
of the Dirac mass matrix ([T9|). As a result it is given by 



y = (g + + v t) + l/e(z/ M - u T ) 

a >/(* + (a: + y)e) 2 + (z + (s-y)e) 2 

The state orthogonal to it is 

v = (z + xe)(u^ - u T ) - yeju^ + v T ) 
" ^{z + {x + y)ef + {z + {x-y)ef 

We concentrate on the case that the dominant contribution to the second heaviest state is due 
to the exchange of the right handed neutrino with mass rri2 which couples to the second row in 
the Dirac mass matrix. The mass eigenstates are given by 



^3 ^ Va (24) 



\z\e li u e + \r\u b 

V2 ^ . = ; = (25) 



\r\p e — \z\e ^Uf, 



where 



\J z 1 +r 2 



^-{zu - a u wy) (27) 



z 



Note the importance of the higher order terms in determining the lighter neutrino eigenstates. 
In particular, the SU(3) family symmetry aligns the leading terms in the (2, 2) and the (2, 3) 
directions and the (1,2) and the (1,3) directions; at this order b' u = c' u = z, d' u = a' u = aw so 
the coefficient a' y b' y — d v d! v vanishes in leading order (this is the second case mentioned above). 
For this reason the second heaviest neutrino can have a large v e component due to the higher 
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order terms. Let us consider the phenomenological implications of this fit. The mixing angle #23 
relevant to atmospheric neutrino mixing is given by 



tan 2 6> 



2: J , 



z + (x + y)e 



z + (x — y)t 



(28) 



Note that it is relatively insensitive to the higher order terms of eq(|T§D. 
The mixing angle 9u relevant to solar neutrino mixing is given by 



tan 2 # 
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(29) 



With these results we may now ask whether it is possible to obtain near bi-maximal mixing with 
the parameterisation of eq([T9"D, with the parameters constrained to fit the quark masses. As the 
latter do not determine all the parameters precisely, we will use the remaining freedom to try to 
generate all the phenomenologically allowed cases, namely the LMA, LOW and VAC solutions 
although we note that the recent SNO data prefers the LMA solution. 



5.2.1 Solution for k=0 



z 


0.91 ±0.06 


X 


-(2.95 ±0.1) 


y 


0.74 ±0.07 


4> x 


0.05 ±0.01 


<Py 


0.19 ±0.01 


W 


0.57 ±0.05 


U 


-(0.55 ±0.12) 


e 


0.21 ±0.01 


e 


0.07 ±0.01 


t 2 


0.58 ±0.02 


t 2 
1 23 


1.49 ±0.05 


(mA 2 < 

\M 3 J ^ 


(0.8 ±0.4) x 10~ 5 


a u 


2.4 ±0.3 



Table 1: The parameters for a fit to the down quark mass matrix under the LOW constraint for 
9 12 in Table 6. 



For k — 0, i.e. a v = —3, the condition assumed above that the right handed neutrinos 
of mass mi and mi dominate the see-saw contribution to the heaviest and next heaviest light 
neutrino masses respectively is > 6 . In this case the masses of the light neutrino mass 

L J mi 7T12 

eigenstates are given approximately by 

M 3 = — 2z 2 v 2 (30) 
mi 

e 6 

M 2 = — z 2 (l + tan~ 2 9 l2 )v 2 (31) 



m 2 

Mi < — v 2 (32) 
m 3 
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z 


0.95 ±0.01 


X 


-(3.15 ±0.31) 


y 


0.64 ±0.03 




0.039 ±0.006 




0.19 ±0.01 


W 


0.53 ±0.03 


u 


Id K 4- n 1 "\ 
— ^U.O ± O.L) 




22 ± 01 


e 


0.07 ±0.01 


t 2 


1.27 ±0.04 


t 2 
L 23 


1.4 ±0.11 






(MzY < 

\M 3 J 


(6.0 ± 1.6) x 10~ 6 


a u 


2.7 ±0.1 



Table 2: The parameters for a fit to the down quark mass matrix under the VAC constraint for 
0io in Table 6. 



Thus for — > we have 

mi ni2 

As we have discussed a large value for 9\2 appears naturally. However M2/M3 is of 0(e 2 ) and for 
choices of the parameters z, uj of 0(1) and #12 in the observed range, we can not obtain a large 
enough value for the ratio M 2 /M 3 to fit the LMA solution. Thus in this case we can only obtain 
the LOW and VAC solutions of the data analysis of neutrino oscillations and mixings. The 
remaining parameters are constrained by the requirement that the parameterisation of equation 
fll9] ) fits the down quark mass matrix (the first fit of eq(f|)). The results are presented in Tables 
H] and H]. These may be compared to the fits to neutrino masses and mixing angles coming from 
the analysis of neutrino oscillations; for convenience we summarise the present situation in the 
Appendix in Tables || ^,^| and [l(| As may be seen we can obtain a good description of both the 



M2 f z Y ( l ±tan" 2 12 \ 2 



LOW and VAC solutions. 

An important cross check of the structure of the GUT relations used in this analysis comes 
from the up quark mass matrix. It is given by eq(|i~9D with e = e and a = a u = 1 (since k — 0). 
The prediction a u = 1 is the analogue of the prediction b' = b in eq@. Given the parameters z 
and uj from the fit of Tables [l] and || we can test this prediction by using the up quark masses 
to determine a u . The results are also given in the Tables. We may see that the inclusion of the 
higher order corrections has actually made the fit to the up quark mass matrix worse than the 
case (b' = b) tested in eq(P|) casting doubt on the viability of the k — solution. 



5.2.2 Solution for k=2 

For k = 2, i.e. a v = 0, the condition that the right handed neutrinos of mass mi and ni2 
respectively dominate the see-saw contribution to the heaviest and next heaviest light neutrino 
eigenstates masses is — > O(-). In this case the masses of neutrino mass eigenstates are given 

mi v mo ' 00 
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z 


U. (0 ± U.14 


X 


( 1 op _i_ n ci \ 

— (l.ob ± U.51J 


y 


1 97 _i_ n 1 


i 

<Px 


0.14 ± 0.04 


<Py 


0.19 ±0.05 


W 


0.77 ±0.06 


U 


-(0.87 ±0.15) 


€ 


0.18 ±0.02 


e 


0.06 ±0.01 


t 2 


0.4 ±0.1 


t 2 
L 23 


1.7 ±0.2 






(Mz\ 2 

\M 3 J 


(3 ±0.H) (^) 2 


a u 


1.6 ±0.4 



Table 3: The parameters for a fit to the down quark mass matrix under the LMA constraint for 
012 in Table 6. 



0.95 ±0.05 
-2.80 ±0.13 
1.37 ±0.16 
0.09 ±0.01 
0.18 ±0.01 
0.77 ±0.02 
-0.86 ±0.06 
0.58 ±0.04 
1.75 ±0.18 
0.19 ±0.03 
0.07 ±0.01 

1.77 ±0.11) 
1.9 ±0.1 



Table 4: The parameters for a fit to the down quark mass matrix for the solution k = 2, under the 
LOW constraint for #12 • 
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z 


0.96 ±0.06 


X 


-3.20 ±0.13 


y 


0.87 ±0.13 




0.05 ±0.01 


<f)y 


0.19 ±0.01 


in 


61 ± 02 


u 


-0.61 ±0.04 


t 2 


1.23 ±0.01 


t 2 


1.55 ±0.20 


e 


0.20 ±0.03 


e 


0.07 ±0.01 






( M 2 \ 2 
\M 3 J 


(0.82 ±0.1) (^) 2 


a u 


2.4 ±0.2 



Table 5: The parameters for a fit to the down quark mass matrix for the solution k = 2 and under the 
VAC constraint for 9%2 ■ 



approximately by 



e 6 

M 3 = — 2z 2 v 2 (34) 
mi 



€ 



6 



M 2 = — (2u 2 + z 2 W (35) 
m 2 

Mi < —w 2 (36) 
m 3 



M 2 ^ 2m 2 + z 2 (m x 
M3 2^ \m 2 

Since the condition > O(^) is readily satisfied for mi/m 2 < 0(1) the LMA solution for 
neutrinos can readily be reproduced through a choice of the ratio mi/m 2 . The results of the fit 
to the down quark mass matrix is given in Table [| This solution is obtained when tan# 2 2 = 0.4, 
which is the favoured value for the LMA solution. As may be seen from a comparison with 
Tables [| [F] and ||, the agreement with neutrino mixing angles is good. 

For the case k = 2 the prediction for the up quark mass matrix is that a u = 2. We may see 
from Table || that this is in better agreement with the up quark masses than the case k = and 
also with the case tested in tested in eq([|) when the effects of higher order operators were not 
included. This is an encouraging indication that the case k = 2 may be realised. 

It is also possible to obtain the LOW and VAC solutions for k = 2. The results are presented 
in Tables f| and [5]. These may be compared to the fits to neutrino masses and mixing angles 
coming from the analysis of neutrino oscillations in good agreement with the fits to experiment 
which are summarised in the Appendix in Tables 1710 and 110. In contrast to the k = case, 
the value for a u is in excellent agreement with the measured up quark masses. 
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5.3 The prediction for Vi3 

Our analyses has so far investigated the implications for V23 and V12 which follow from an 
enhanced flavour and family symmetry. What are the implications for V13? This is a particularly 
interesting question for only if V13 is quite large will there be any prospect of seeing CP violation 
in future long baseline neutrino experiments. Due to the form of the hierarchy in the mass 
matrices, both the charged leptonic mass matrix and the effective neutrino mass matrix, it is 
possible to diagonalise them by the rotations R23R13R12 together with diagonal matrices carrying 
the phases. From this the 13 (e3) V MNS mixing matrix is given by [13, |TJ 



rrMNS _ v I I iui _ I J v I iuJ 2 i J „Ll v iuj 3 (no\ 
v 13 ~ i 13 c 12 c 13 e 6 13 c 12 c 13 c 23 e ^ 6 12 i 23 c 13 e \°°J 

where s^, for i, j = 1, 2, 3 and f — l,u represent the sines of the mixing angles, analogously for 
the cosines and the phases i = 1,2,3 are functions of the phases appearing in the mass matrix. 
The angle S23 is given by |c2 3 s 2 3 — S2 3 c 23 e laJ4 |. Since the neutrino oscillations experiments are 
not sensitive to Majorana CP phases, the standard parametrisation used for the case of quarks 
16 j is commonly used for the case of leptons. In this parameterisation V13 = srse~ tS , where 
S13 should be identified with the sine of the CHOOZ angle and S is the analogue to the quark 
CP violation phase. Thus identifying s 13 with the absolute value of eq(p8[) we can trace the 



contributions from charged leptons and neutrinos to the CHOOZ angle. 

We can see that, in both the k = and k = 2 cases discussed above, the value of s± 3 is 
negligible because of the hierarchy of the Majorana masses we have taken (1713 » 1712, mi). The 
mixings in the leptonic sector are small and are given by 



Wt 1 Vi 1 23 I VI I _Li2_i23 

J ~ i33_. J ^33 /OQ\ 



|r 22 Y 1 — 1 ''' 

^33 



I 
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where Y\ 2 are the elements of the charged lepton mass matrix. From eq (|l9D we see that S 
0(e 3 ) and s[ 2 = 0(e 2 ), therefore the latter is the dominant contribution. This result is interesting 
because Yj 2 /Y 2 l 2 ps \pmjm~ tl and hence from eqs(^) and (|39[) 



which is a testable prediction. This gives a value V^ NS ~ 0.07, close to the present bound and 
large enough for the future CP violating experiments to be viable. This point has been made 
explicitly by King [|lf| and the value given is consistent with the results of |17||18l. When the 



third term in equation (|38| ) is dominant, 5, the analogue to the quark CP violation phase, is 
— co> 3 , which in terms of the elements, \Y? |e 1 ^, of the effective neutrino mass matrix it is given 
by: 

5 = - U3 « 7 ^ 3 -I12 zJk (41) 



where 



^ ~ F^.cos^-^+^cos^-^ 1 } 

723 ~ 023-033 (43) 
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for ij = 12, 13 and k ^ In the same way as discussed for the quarks, Section |2|, we may 
choose the phases to be <p l , the phase of the (1, 2) element of the Dirac M v and x l the phase of 
the (1, 3) element of M l , then 5 ~ <p l . 

The process of extracting the angles and phases from successive rotations is general and can 
be applied directly to hierarchical mass matrices; for a further work on this applied to leptons 



and the implications for CP violation see pA 
5.3.1 n — > cy 



Several groups [TT3, have pointed out that off diagonal lepton Yukawa couplings can lead to 
unacceptably large contributions for lepton family number violation processes such as fi — » &y. 
In a supersymmetric framework the branching ratio, BR{^i — > e^y) has the form 

BR(n -»• erf) <x | (yj In ^Y v ) 21 1 2 tan 2 /?, (44) 

where Y^, is the matrix of Yukawa couplings for Dirac neutrinos (in the basis in which charged 
leptons are flavour diagonal), Mx is the scale of Grand Unification, Mr is the scale where right 
handed neutrinos decouple and tan /3 is equal to the ratio of the vacuum expectation values of the 
two Higgs doublets of the MSSM. Since BR^fi — > ej) depends as well on the soft supersymmetric 
mass spectrum is useful to determine instead the matrix element C^ e = (Yj \n(M x / M R )Y v ) 2 i 
and hence study the predictions of the structure of Y v and the scale of Mr for lepton flavour 
violating processes ||20| . 



We can compute the element C Me for our texture of Y v . For the preferred case k = 2 we have 
C Me w 6 x 10~ 3 for the LMA solution, C^ e w 5.6 x 1(T 3 for the LOW solution and w 4.2 x 10" 3 
for VAC solution . Even for the case that tan/3 is large (h t ~ hf,), due to the hierarchical form 
of the Yukawa couplings, eq(|19D, |(Fjl^)2i| tan f3 is small and therefore does not conflict the 
bounds for jj, — > ej f\. As can be seen from Figure 2 [2(J the coefficients C^e given above fall 



below current experimental bounds lyj for a wide range of the soft supersymmetric breaking 
parameter space 0. 

5.4 Summary and conclusions 

Due to the see-saw mechanism the significant differences between quark and lepton mixing 
can be explained while keeping the form of their Dirac mass matrices the same. As a result 
the observed masses and mixings can be accommodated in a theory in which there is a very 
large underlying family and flavour symmetry group. We have explored the phenomenological 
implications of an SU(3) family symmetry together with a 5*0(10) GUT flavour symmetry and 
additional Abelian family symmetries, chosen to restrict the allowed Yukawa couplings. Allowing 
for spontaneous (perturbative) breaking of this group we found a symmetry breaking scheme in 
which the observed hierarchical quark masses and mixings are quantitatively described together 



7 The structure of the effective mass matrix for low energy neutrinos (provided by the see-saw mechanism) in 



the models considered here belongs to the Class 3 of reference [20 . 

8 To compare with the bounds presented in |20[ we need to multiply the upper bounds for C^ e by 10/50, since 
tan/3w50. 
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with the hierarchy of charged lepton masses and an hierarchical structure for neutrino masses 
with near bi-maximal mixing in the lepton sector. In this the presently unknown mixing angle, 
(#M7vs)i3; is determined mainly by the mixing in the lepton sector. While smaller than the other 
mixing angles, it is close to the present limits, and is large enough to allow for significant CP 
violating effects to be visible in future long baseline neutrino experiments. 

Given the very large underlying symmetry, the fermion masses are heavily constrained. The 
perturbative breaking ensures an hierarchical structure for the masses and in terms of the break- 
ing parameters, the order of magnitude of the ratios of the quark and lepton masses is determined 
once the family symmetry properties of the fields are determined. For quarks, having constrained 
the family charges to fix the down quark mass ratios, one obtains one order of magnitude pre- 
diction for the ratio of up quark masses which is in good agreement with experiment. Further 
the family symmetry ensures an approximate texture zero in the (1,1) matrix element which 
gives rise to the successful GST relation between the CKM mixing and the mass ratios of the 
first two generations of up and down quarks. 

The extension of the family symmetry to leptons can be done together with via an enlarged 
GUT flavour symmetry, 5*0(10). This results in predictions for the charged lepton masses in 
terms of the symmetry breaking pattern of 50(10). For one particular choice, the predictions are 
in good agreement with measurement, reproducing the Georgi Jarlskog structure for the light 
lepton masses. In addition the 50(10) symmetry relates the Yukawa couplings of the up and 
down quark sectors and replaces the order of magnitude prediction for the ratio of up quarks 
by an absolute prediction. For the preferred choice of symmetry breaking this prediction is in 
excellent agreement with experiment. 

The neutrino mass matrices are also strongly constrained by the symmetry. If the expansion 
parameter for the Majorana masses is smaller than that for the up quarks, the Majorana mass 
matrix is approximately diagonal in the family symmetry basis. In this case the lepton mixing 
angles are determined by the Dirac mass matrices for the charged leptons and neutrinos and 
these in turn are related by the GUT symmetry to the Dirac mass matrices of the quarks. 
As a result the lepton mixing angles^] are determined by the properties of the quarks. The 
mixing angle relevant to atmospheric neutrino mixing is well determined by the leading order 
operators in the symmetry breaking expansion. Through vacuum alignment in the 577(3) family 
sector this mixing angle is large and consistent with the measure atmospheric mixing angle. 
The solar mixing angle turns out to be quite sensitive to subdominant terms in the symmetry 
breaking expansion. As a result one obtains only an order of magnitude prediction for this 
angle determining it to be of 0(1). For the favoured 5O(10) symmetry breaking pattern the 
coefficients of the sub-dominant operators needed to obtain a quantitative agreement of this 
mixing angle with data are consistent with the constraints on the expansion which follow from 
fitting the down quark mass matrix structure. Finally the value of (0mns)i3 is determined to 
be approximately ^Jrrie/m^, coming mainly from the charged lepton sector. 

To summarise, the data on all quark and lepton masses and mixings is qualitatively consistent 
with a significantly enlarged family and flavour symmetry. Given the disparity between these 
quantities in the quark and lepton sectors this is already remarkable. On a more quantitative 
level, a specific pattern of family symmetries and symmetry breaking leads to a quantitative 
prediction for seven of the parameters of the Standard Model, namely {Qckm)v2, nT'u' m t/' m ci 

9 The neutrino masses, while constrained in some case, are largely determined by the unknown Majorana mass 
eigenvalues. 



19 



m e , m M , m T and (9mns)23 together with an order of magnitude prediction for (9mns)i3- Such 
agreement is encouraging and suggests one should try to construct a complete underlying theory, 
be it a SUSY GUT or perhaps a superstring theory. 
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A Appendix 





LMA 


LOW 


VAC 


Ref. 


tan 2 6i2 

\M 2 \ 2 
\Ms\ 

(tan 2 6>i 2 )ex P 

\Am 2 12 \ 

l Am !3lcxp 
\Am 2 12 \ 

l Am 2 3 lexp 


(4.0 ±1) x lO" 1 
3^ 

m 2 

4.0 x lO" 1 
(0.3,33.6) x 10~ 2 

(0.5,30.8) x 10" 2 


(5.8 ±4) x lO" 1 

2 

mi 

5.8 x 10" 1 
(0.48,10) x 10" 5 

(0.73,10) x 10" 5 


(1.23 ±0.1) 
8.2 x lO" 1 ^! 

(0.68,1.8) 
« 1.16 x 10~ 7 

« 1.44 x 10" 7 


[22, H, 24| 
25] 

HI 



Table 6: Comparision between the predictions for k = 2 and the experiments for tan 2 # 12 and 
| -||- 1 2 for the LMA, LOW and VAC solutions for solar neutrinos. The confidence levels presented 
are at 3a. 



21 



Atmospheric+chooz information 


Am 2 2 




t 2 

L 23 




t 2 




Ref. 


BFP 


Range 


BFP 


Range 


BFP 


at 90 %CL 




3.1 x 10- 3 
2.5 x 1(T 3 


(1.1,7.3) x 10" 3 

at 99 %CL 
(1.2,4.8) x 10~ 3 
at 3cr 


1.4 
1 


(0.33,3.7) 
at 99 %CL 
(0.43,2.3) 
at 3a 


0.005 


< 0.055 


|25| 
|26| 



Table 7: The constraints on neutrino mixing parameters coming from Atmospheric neutrino data 
and from CHOOZ. The second analysis cited has been performed with updated atmospheric data. 
For different values of tf 3 there are different two dimensional CL. regions for the variables £23 
and Arn| 3 ; here we present the appropriate for the value of t 2 3 obtained with the mass matrix 
structure discussed here. 



LMA, solar information 


Am{ 2 (eV) 2 




t 2 




g.o.f. 


Ref. 


BFP 


at 3(7 CL 


BFP 


at 3(x CL 






5.5 x 10~ 5 

5.6 x 10~ 5 
6.2 x 10~ 5 


(0.23,3.7) x 10" 4 
(0.22,2.2) x 10~ 4 
(0.23,3.4) x 10~ 5 


4.2 x 10- 1 
3.9 x 10- 1 
4.0 x 10- 1 


(2.4,8.9) x 10- 1 
(2.0,6.4) x 10- 1 
(2.3,7.9) x 10- 1 


49% 
84% 




23 
27 




4.5 x 10~ 5 


at la CL 
(3.1,7.2) x 10- 1 


4.0 x 10- 1 


at la CL 
(3.2,4.8) x 10" 1 
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Table 8: Experimental constraints for LMA solution. BFP refers to the best fit point given in 
the cited references, the 3cr=99.73% CL or the la CL. is shown. 



LOW, solar information 


Am 2 12 


Am 2 12 


t 2 


t 2 


g.o.f. 


Ref. 


BFP 


at 3a CL 


BFP 


at 3<7 CL 






1.0 x 10~ 7 

1.1 x 10~ 7 


(0.5,1.1) x 10- 7 
(0.6,1.2) x 10~ 7 
(0.9,1) x 10~ 7 


7.1 x 10- 1 
6.9 x 10- 1 


(5.5,10) x 10- 1 
(4,9) x KT 1 
(6,7.5) x 10- 1 


45% 
69% 


m 
m 



Table 9: Experimental constraints for the LOW solution. BFP refers to the best fit point given in the 
cited references, the 3<r CL has been estimated from those references. The analyses cited have been 
performed using the latest SNO information. 
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VAC 


Am 2 12 


Amf 2 


t 2 


t' 2 


g.o.f. 


Ref. 


BFP 


at 3cr CL 


BFP 


at 3a CL 






4.6 x 1(T 10 


(3.5,5.7) x 10" 1U 


2.4 x 10° 


(0.3,3.5) 


42% 


|2S| 



Table 10: Experimental constraints for VAC solution. BFP refers to the best fit point given in the 
cited references, the 3a CL has been estimated from those references. The analysis cited has been 
performed using the latest SNO information and with an enhanced CC cross section for deuterium, as 
quoted in pS|]. 
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